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We consider systems with slab geometry of finite thickness L that undergo second order phase 
transitions in the bulk limit and belong to the universality class of 0(n)-symmetric systems with 
short-range interactions. In these systems the critical fluctuations at the bulk critical temperature 
0^ , Tc,oo induce a long-range effective force called the "thermodynamic Casimir force". We describe 

' the systems in the framework of the 0(n)-symmetric ^''-model, restricting us to the large-n limit 

n oo. In this limit the physically relevant case of three space dimensions d — 3 can be treated 
analytically in systems with translational symmetry as, e.g., in the bulk or slabs with periodic 
5— ( . or antiperiodic boundary conditions. We consider Dirichlet and open boundary conditions at the 

surfaces that break the translational invariance along the axis perpendicular to the slab. From 
the broken translational invariance we conclude the necessity to solve the systems numerically. We 
evaluate the Casimir amplitudes for Dirichlet and open boundary conditions on both surfaces and 
for Dirichlet on one and open on the other surface. Belonging to the same surface universality class 
we find the expected asymptotic equivalence of Dirichlet and open boundary conditions. To test 
the quality of our method we confirm the analytical results for periodic and antiperiodic boundary 
' conditions. 

o 

D ' PACS numbers: 05.70.Jk, 68.35.Rh, ll.lO.Hi, 68.15.-l-e, 75.40.-s 

^ ^ Keywords: Casimir effect, fluctuation-induced forces, scaling functions, large n-limit, spherical model 
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I. INTRODUCTION 



. In a fluctuating medium confined between two surfaces long-range forces may arise. This effect was first predicted by 

■ H. B. G. Casimir [1] who considered a confinement of the electromagnetic vacuum fluctuations between two perfectly 
I , conducting metal plates. The boundary conditions imposed by the two plates modify these fluctuations in such a 

■ way that the energy of the system becomes dependent on the separation of the plates. From this energy gradient in 
C ' the separation a long-range force, the so called Casimir force arises. This force was experimentally measured for the 

flrst time in 1998 by Mohideen and Roy [2]. After this breakthrough more measurements followed [3, 4], conflrming 
Casimir's prediction with high accuracy. 

An analogous effect is the so called thermodynamic Casimir effect which was first proposed by Fisher and de Gennes 
' [5]. It is caused by critical fluctuations of a medium near its bulk critical point Tc^oo, where these fluctuations are 
^ , correlated over an infinite range. The spatial confinement of such fluctuations also gives rise to a long-range Casimir 
' force. 

\^ , The thermodynamic Casimir effect was experimentally proven for the first time by Garcia and Chan in thin films of 
' liquid "^He at its transition temperature to superfluidity [6] . They measured a characteristical dip of the fllm thickness 
■ ' _ at the critical temperature, that can be attributed to the thermodynamic Casimir force. 

. The experimental data are excellently reproduced by Monte Carlo simulations of the XY-model with open boundary 
' conditions [7] which is in the same universality class as the superfluid transition in ^He. 

: Due to the presence of strong fluctuations in critical systems there are only a few exact results for the thermodynamic 
^ Casimir effect yet. These are usually restricted to the case of periodic boundary conditions along the slab. Though 
, there are perturbation theory results up to second order for experimentally relevant boundary conditions as, e.g., 
Dirichlet boundary conditions [8], exact results for these are still lacking. 
5h The remainder of this paper is organized as follows. In the next section we carry the elementary features of the 

. - - 1 thermodynamic Casimir effect together, and introduce the large-n limit. In Section II the model considered in the 
present work will be defined and its excess free energy will be derived in Section III. This calculation leads to a 
self-consistent eigenvalue problem, which is solved numerically, providing the possibility for numerical estimation of 
the Casimir amplitudes. The paper is closed with a discussion of the achieved results. 



A. Casimir effect 



As predicted by finite-size scaling theory a statistical system confined between two parallel surfaces exhibits a 
dependence of the free energy on the separation L of the surfaces [9] . This L-dependence leads to an effective force 
between the surfaces which is a typical example of a finite-size effect. This thermodynamic Casimir force is defined 
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as [5, 10-12] 



rc{T,L) = -kBT^^^^^, (1.1) 



where /ex denotes the reduced excess free energy per unit area, which is given by /ex = /l — Lf\^ with the total free 
energy per unit area //, and the reduced bulk free energy density /b. The limit A ^ oo oi the cross sectional area 
has been taken. 

For temperatures away from the bulk critical point the bulk correlation length £,oo{T) is much smaller than the 
macroscopic separation L. In this regime the Casimir force decays exponentially for large L. 

As it is generally known, ^oo(T) diverges as jT— Tc.ool at the bulk critical temperature Tc.oo, with critical exponent 
V, so in the vicinity of T^^oo the bulk correlation length is of the same order as L. This implies that the Casimir 
force Tc{T, L) extends to distances much larger than the microscopic scale a (~ lattice constant). In this regime the 
Casimir force decays algebraically for large L [5] . 

Examining the excess free energy more precisely one finds that it can be decomposed into two contributions 

/ex(T, L) = /,(T) + /res(T, L) . (1.2) 

The first one is the _L-indcpendcnt nonlocal surface excess free energy fs{T) = fcx{T, oo) [13]. The second one is the 
residual finite-size contribution fj:cs{T,L) which contains the whole information about the L-dependencc of the free 
energy. 

At the bulk critical point Tc,oo the finite-size contribution has the characteristic algebraic behavior [31] 

/res(Tc,oo,i) ~ Ac L-^'^-^\ (1.3) 

L—*oo 



leading to a long-range effective Casimir force 



{d-l)AcL-'^. (1.4) 



The quantity is an universal finite-size quantity, called the Casimir amplitude [5] . Universal means that it depends 
only on the bulk universality class of the phase transition and the surface universality class of the boundary planes 
[13]. It is in particular independent of microscopic details of the system, such as, e.g., the lattice structure. 
Near the bulk critical point the residual free energy density and the Casimir force obey the scaling forms 

/,es(r, L) ~ e{L/U{T)) (1.5) 

and 

^^^I^^L-'^^iL/^^iT)), (1.6) 

where Q{x) and ^{x) are universal functions with the scaling argument x = L/^oo- At Tcoo, where L/^oo = 0, they 
take the values 6(0) = Ac and S(0) = (d - l)Ac. 

There is no complete analytic theory for the thermodynamic Casimir effect joining the regimes above and below 
Tc^oo- In the framework of the 0'*-theory perturbative calculations up to second order for different boundary conditions 
at and above Tc^ao [8, 14-16] were performed. There are also exact results for the spherical model with periodic [12] 
and antiperiodic [17] boundary conditions. However, for slab systems with non-periodic boundary conditions no exact 
results are available. This is the aim of the present work. 



B. The large-n limit and the spherical model 



To investigate the thermodynamic Casimir effect one needs a model that shows the characteristic critical properties 
of a second order phase transition. There are various classes of models that exhibit critical behavior but here we 
restrict ourselves to the class of classical 0(n)-invariant n-vector lattice models. Those can be represented by the 
Hamiltonian 

<i,j> 
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where the S = (si,...,s„) arc n-component vector spin variables, restricted to the local constraint \S\'^ = 1. The 
exchange coupling J is assumed to have a positive nonzero value and the sum runs over all nearest neighbors. 

In the limit n qg, which is called the large-n limit this model can be solved analytically for the bulk or for 
systems of finite size [18]. This case is connected to the so called the spherical model, which was first introduced and 
solved by Berlin and Kac [19] and is defined by the Hamiltonian 

<i,3> 

Here the Si are scalar spin variables fulfilling a global constraint ^isf = N called "spherical constraint" , where N is 

the number of lattice sites. 

This model is even more easy to handle if one uses the constraint (J2i ^i) = ^ where one takes the statistical mean 
over the sum. It is then called the "mean spherical model" but the both formulations are equivalent in translationally 

invariant systems. 

It was found by Stanley [20] that these spherical models are equivalent to the large-n limit for translationally 
invariant systems as e.g. the bulk case or systems with periodic boundary conditions. 

To hold this equivalence in non-translationally invariant systems the spherical constraint has to be modified [21]. In 
the case of a slab where the boundary conditions are non periodic one needs a spherical constraint for each layer parallel 
to the surfaces. These multiple constraints make the analytical treatment impossible so that numerical calculations 
are needed. This is the challenge we are going to face here. 

The large-n limit can also be understand as the zeroth order of a systematical expansion in 1/n [22] . This connects 
our theory to the physical relevant case of e.g. three order parameter components known as the Heisenberg model. 
But this expansion is more useful for formal studies than to comparison to experimental observations. 



II. DEFINITION OF THE MODEL 



In the last chapter we have discussed 0(n)-symmetrical lattice models but now we turn to the description of the 
universal critical behavior in the framework of an 0(n)-symmetrical (^'^-model. The connection between both is that 
the n-vector lattice model can be mapped on the n- vector (^^-model by means of coarse graining [23]. 

We consider a d-dimensional slab of finite thickness L, which occupies the volume 2J = x [0,L]. Let the Xj 
with j = 1 . . . , d, be cartesian coordinates and we write z to denote the coordinate of the finite direction called 
the perpendicular direction. The full position vector is given by x = {y,z), where y = {xi, . . . ,Xd-i) denotes the 
coordinates of the directions along the slab called parallel directions. 

As discussed in detail in [13] the general Hamiltonian of 0^-models with surfaces is given by 

n[<P{x)]= f C^[<P{x)]dV+ [ Ca^[<p{x)]dA. (2.1) 

The densities C<x}{x) and Ca<xs{x) depend on the n-component order parameter field 0(a;) and its derivatives. This 
order parameter field is a classical n-component vector field that, using the language of magnetism, can be interpreted 
as the local spin density. The first part is the bulk Hamiltonian and the second is the additional surface Hamiltonian. 
The latter is taken in account because the coupling at the surface might deviate from that in the bulk. 

The bulk density £«o(a;) is given by 

a— 1 

where f is the temperature variable often called mass and u is the coupling constant. 

Note that in contrast to (/(^-theories with a finite number of order parameter field components there is an additional 
1/n in front of the (p'^-term, which is needed to make the limit n — > oo well defined [24] . 

The general boundary density £dm{x) reads 

£ommx)] = ^<p'{x). (2.3) 

The c{x) is the sm-facc enhancement that controls the additional couplings at the surface. It may have different values 
on the two boundaries i9*8i and 0^2, i-e., 

' forxe9»i, 
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In our framework these boundary terms need not to be taken into account any longer because they give no contri- 
bution to the finite size energy that we are interested in and one can take account for the boundary conditions on the 
surface without those additional terms. The boundary conditions will be taken into account in the calculation of the 
eigensystem of the system. How this is done will be described in the next section. 



III. THE FREE ENERGY 



Wo now want to derive the exact expression for the free energy of the system. In the large-n limit such an expression 
can be obtained in different ways. One can use e.g. variational techniques as well as the saddle-point approximation 
[25]. Both become exact in the limit n ^ oo. In the following we resort on the saddle-point approximation method. 

We start the calculation by computing the canonical partition function Z of the system. It is given by the functional 
integral over all possible configurations of the order parameter field 4*{x). We write 

Z = I P[0]e-«[<^1 = I V\4>\ exp |- y A'^x \ j^iVM' + \<i>" + ^ } ' (^-l) 

To determine the partition function with the saddle-point approximation we transform the last expression by 
applying the Hubbard Stratonovich decoupling technique [26] and introducing an additional parameter ^j). This leads 
to 



Z = C j V[4>\ J P[zV']exp|-i J d'^x 



(3.2) 



The constant pre-factor C is irrelevant since it does not contribute to the excess free energy. The introduced variational 
parameter ip is chosen to depend on z because the boundary conditions break the translational invariance along that 
axis. 

The integration over <p{x) is now Gaussian and can be carried out (see e.g. [27]). This leads us to the expression 



Z = C 



exp 



-|TrIog(-V2 + f-^(^))-gA / dzV'(^) 



(3.3) 



The remaining functional integral can be evaluated by the saddle-point approximation. This means that we take 
the value of the integrand at its saddle point as solution of the integral. The number of order parameter components 

n serves as large parameter in the exponential function that makes the saddle-point a sharp peak as required for this 
approximation. With growing n the peak becomes sharper so that the approximation gets better and will become 
exact in the limit n — > oo. So we conclude that the partition function is given by 



Z = C exp 



3n 
2^' 



Jo 



(3.4) 



Prom the partition function we obtain the reduced free energy using its definition F{L) = T /k^T = — \ogZ. It is 
useful to consider the reduced free energy density per cross sectional area and order parameter component, /l = F/An 
with oo, that comes out as 



/l = /o + 2 



(d-i) 



Tr, \og{-dl + p2 + f - V'o(^)) - ^ / ^l^^) ■ 



(3.5) 



Note that we have set ■i/' = V'o in equations (3.4) and (3.4) where is the saddle-point value. It is defined by 
requiring that at the saddle-point the integrand must be extremal in i/', and obeys the self-consistent equation 



V'o(2) 



[d-i) 



Tr, 



1 



(3.6) 



We now examine this self-consistent equation and determine the saddle-point value and the eigensystem along the z 
axis. 

But first we want to remark the following: The problem can also be approached with the help of diagrammatic 
perturbation theory. From this theory one obtains the Dyson equation 



G(fc) =k + T-t{k) 



(3.7) 
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for the two point correlation function of the system often named propagator G{k). The contributions from aU orders 
of perturbation theory are summed in the mass-operator S(fc). UsuaUy this cannot be done without approximation 
but, in the large-n hmit the contributions to the mass operator can be summed up exactly (see [24, 28]). That 
mass-operator has the same self-consistent structure as ipo{z) in equation (3.6). So in comparison we can identify 
^'0(2) = S(-2) and we will use S(^) instead of V'o(-z) the following. 

For simplicity we write the mass-operator T,{z) in the mixed p-2:-representation 

n^)=-\ rv'^^^^p!- (3-8) 



Here we define 

r-(d-l) 



dp. 



i=l 



=0 2r 



and use this notation in the following. 

The <fv{z) and the Cj/ are the eigenfunctions and the eigenvalues used to expand the ^-dependence and are given 
by the equation 

[-92 + S(^)] ^^z) = e.^u{z) . (3.10) 

Obviously the quantity Ti{z) is not well defined because the p-integration is UV-divergent. Usually one would have 
to deal with the procedures of regularization and rcnormalization to give the divergent expression a physical meaning. 
But our case is much simpler because we are not interested in the mass-operator itself. Thus we define another 
quantity which is UV-convergent by taking the difference of ^{z) and a quantity which is UV-divergent to the same 
order. 

For this purpose we decompose the mass f into f = fc + Sf. The fc gives the value of f at which the corresponding 
bulk system reaches criticality and 6f gives the deviation from this critical value. The fc is given by the expression 
[24] 

n(d-l) 



^r"v°°^^— (3 11) 



It is UV-divergent to the same order as ^{z). From the sum of both we now obtain a UV-convergent expression. 
We define the new quantity 

Uiz) = E(.) + fc = ^ ( V ^-(-)f(-l - r . (3.12) 

Introducing this quantity the eigensystem is given by 

[-dl + U{z)\ Mz) = ^I'Mz) ■ (3.13) 

To simplify the further discussion we will interpret the two terms on the right hand side of equation (3.12) analogue 
to Rcf. [28]. The first term on the right hand side is the propagator of the slab system for temperature deviations from 
the critical point defined by Sf. The second term is the propagator of the corresponding bulk system at criticality 
where Sf = 0. With this definitions one may write equation (3.12) in the following form 

Uiz) = -J (Gip; z, z) - Gbuik(p; q)) (3.14) 

Equations (3.12) and (3.13) form a self-consistent eigenvalue problem similar to the one dimensional Schrodinger 
equation of a particle in a potential U{z). The numerical solution of this problem is the content of the following 
subsection. 
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FIG. 1: Numerical data of the potential for a system with Dirichlet boundary conditions. 



A. Discretization of the problem 



Using a numerical method we now calculate the potential Vi^z) and the corresponding eigensystem self-consistently. 
Therefore the model needs to be discrete along the z-axis. We introduce a discretization of the z-direction with 
points and the lattice constant a in the sense that z — > with zi G {a, 2a, 3a, iVa}. The lattice constant is chosen 
to be a = 1 without loss of generality. 

This discretization leads to some more changes. At first the number of eigenvalues and functions is restricted 
to the number of discrete points iV. At second there is also a change of the dispersion relation in the corresponding 
bulk model. The in the second term of the right hand side of equation (3.12) has to be replaced by the dispersion 
relation 



4 . ^ qa 



(3.15) 



of a lattice theory. Due to the periodicity of (3.15) the integration with respect to q is restricted to the first Brillouin 
zone. 

Implementing these modifications the potential U(z) displayed in (3.12) becomes 



■ N 

E 



Aq 



1 



27rp2 + ^sin^(^) 



(3.16) 



where the corresponding eigenvalue equation (3.13) now reads 

{-dX+U\^, = e,^, (3.17) 

Here the d\ and TJ are finite N x N matrices and the (fi y are the corresponding finite A''-component vectors. We 
specify these objects now: The &], stands for the symmetric differentiation operator for the second derivative of a 
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FIG. 2: Numerical data of the potential for a system with open boundary conditions. 



discrete function. Usually this operator is defined by 

Setting the lattice constant a = 1 explicitly we can represent this operator as the matrix 



— X 


1 


. . 


. 


y 


1 


-2 


1 










1 


-2 













-2 


1 


y 







. 1 


— X 



(3.18) 



(3.19) 



At this stage the boundary conditions of the system come into play. We can control them by choosing adequate 
values for x and y. The reader can convince himself that if one chooses x = 1 [x = 2] and y = this corresponds to 
open [Dirichlet] boundary conditions. If one sets y = 1 [y = — 1] and x = 2 one gets periodic [antiperiodic] boundary 
conditions. 

The matrix U is the diagonal matrix of the U{zi) 

J7 = diag[(7(z,)] (3.20) 
and the eigenfunctions of the discrete model are represented as vectors where the components are given by 

VviZi) = ■ (3-21) 
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FIG. 3: Numerical data of the potential for a system with open boundary conditions on one side and Dirichlet on the other 
side. 



In this discrete form the self-consistent eigenvalue problem can be implemented into an iterative routine to deter- 
mine the eigensystem and the matrix JJ . The calculations were carried out with build-in functions of the program 
Mathematica [29]. 

The iteration starts from an initial guess for the potential. From this guess the associated eigensystem is computed. 
With the eigensystem one computes again the potential for the next iteration step. This procedure is repeated until 
the difference between the steps is small enough to fulfill a certain termination condition. 

Figures 1, 2 and 3 show the numerical results for the potential Uizi) plotted over the scaling variable ZijL. We 
show system sizes L in powers 2" from L — 256a to L = 4096a. 

Figure 1 shows the results for the systems with Dirichlet boundary conditions on both surfaces. The potential as 
anticipated diverges at the boundaries to ensure that the eigenfunctions fall to zero on the surface. In the middle 
of the slab the potential tends with a characteristic algebraic behavior to its bulk value which is zero at criticality 
(see [22]). This tendency is more and more distinct for bigger systems because the influence of the boundaries in the 
middle of the slab decays with increasing system size. 

Figure 2 shows the results for open boundary conditions. In the middle of the slab the potential approaches the 
behavior of systems with Dirichlet boundary conditions. This approach becomes stronger for larger systems. But 
next to the surfaces the evolution differs very strongly from the Dirichlet case. The most remarkable difference is that 
the potential has now a finite, positive value at the surface. 

When there are both boundary conditions mixed in one system as it is shown in Figure 3 we find the diverging 
behavior on the Dirichlet side and the finite values on the open side. The values of the potential in the middle of the 
slab are still converging to the bulk value with increasing system size. 

The Dirichlet and the open system are that similar because they are in the same surface universality class [13]. 
This means that their asymptotical behavior for large L and large values of z far away from the surface is identical. 

To test the quality of our method we have also computed the potential V(^z) for the known cases of periodic and 
antiperiodic boundary conditions with our numerical method. We found that our results are in excellent agreement 
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with the analytical results. A comparison between our results for the Casimir amplitude calculated from this potentials 
and the analytic results is shown in the next section. 



IV. EXCESS FREE ENERGY AND CASIMIR AMPLITUDES 



Now that the U{z) is determined we can calculate the Casimir amplitudes. But first we have to modify the 
expression of the free energy by using the definitions f = fc + (5-f and S(z) = fc — U{z) to have only well defined 
quantities in it. 

As we are not interested in the free energy itself but in the difference of the free energy of the slab and the bulk, 
called excess free energy we need an expression for the free energy of the bulk. We obtain this expression from a 
calculation analogous to that in the last chapter. The translational invariance of the bulk allows one to obtain its free 
energy analytically [22]. 

Now the exact expression for the excess free energy defined in Sec. I can be written in the form 

= I ''l]log(e.+p2 + <5f)-^^''d.(fe-f/(.)f 

r Ad-l) roo 1 or 

- 2/ j ^log(p' + 9' + '5T + [/buik)-^(fc-«7buik)' (4.1) 

Because of the discrete numerical evaluation of the potential U{z) all z dependences in the last expression have to 
be discrete. This is what will be done the next section. 



A. Discrete form of /e 



We now obtain the form of the excess free energy with a discretized z dependence. This means substituting z ^ Zi 
like in Sec. III. With this change we have to replace the dispersion relation of the continuum by that of the lattice 

and restrict the q momentum integrations to the first Brillouin zone (sec also Sec. III). 

Finally the real space integration along the z-axis has to be turned into a sum. This is achieved by following the 



rule/o dz 



From this changes we obtain the excess free energy in the discrete form given by 



/ex — 



(d-1) 



N 



N 



-1) 



P 



6t) 



2u 



.^(fe-[/(za))^ 



qa 



log {p' + ^ sin^(^) + 6t + C/bu 



, 



C^bulk) 



(4.2) 



Since we are interested in the Casimir amplitude we limit the following discussion to the case T 
temperature variable 6f = 0. In the following we use the result that the J/buik equals zero at Tc^oo 
After evaluating the momentum integrals we get the expression 



= Tc^oo and set the 
(see [22]). 



/e: 



1 1 



N 

E 



N 



-,2-d ' 



T^=^U{ia)T{^) 2''-'^^eZ^ , 4a-<^r(f) 



+ 



sec 



~2 



(4.3) 



Examining the term "sec(d7r/2)" in this expression one recognizes that it has a remaining pole term in d = 3 space 
dimensions. This pole term can be imdcrstood by realizing that the g-intcgration is carried out but the corresponding 
summation over the eigenvalues is left unevaluated. The pole term needs to be removed to set d = 3 because the 
evaluation of the remaining sum in arbitrary dimensions is not feasible. 

To make the pole-term visible and isolate it we expand /ex in the space dimension around = 3. In this expansion 
we find a term from order 0{l/{d — 3)) 



N 

E 



2 + ma) 
4a3(d-3)7r 



(4.4) 
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ao 


a2 


£13 


a4 


oo 


0.0297 


-0.0122 


0.0106 


-6.06 


DO 


-0.0466 


-0.0119 


-2.97 


441 


DD 


0.0309 


-0.0115 


-4.53 


427 


PBC 




-0.1530 


-1.23 


51.0 


APBC 




0.2742 


-27.4 


3390 



TABLE I: Fit parameters for the excess free energy, Eq. (4.6). OO means open boundary conditions on both sides, DD means 
the same with Dirichlet boundary conditions, and DO means the mix of both. 



which is singular at = 3. But if we carry out the sums over v an i this pole-term becomes zero in the range of 

numerical accuracy. In other words: The pole-term is just an artifact of the non exact summation of numerically 
calculated objects. The zeroth order of the expansion contains the whole information of the excess free energy at 
d = 3 so we restrict us to the discussion of this term. 

If the lattice constant a is explicitly take to be unity, without loss of generality, we obtain the expression 



1 ^ 

= —Y 



N 



1=1 



+ 



^ N \ - 

e,-J2u{i)] -2 log(47r 

K i=l J . 



N 



e.(l-log(e.))^--^C/^(z) 



(4.5) 



where 7 is the Euler-Gamma constant. In this expression the remaining sums can be evaluated to obtain the excess 
free energies for the given system sizes. 

To compute the Casimir amplitude from the evolution of the numerical evaluated excess free energy with increasing 
system size we fit the data using the ansatz 



/ex = ao + + a^L ^ + a4,L + . 



(4.6) 



where the factor 02 will later be interpreted as the Casimir amplitude. This fitting ansatz is justified by the following; 
From fundamental arguments one knows that the order 0{L^^) is absent and the leading order is the order 0{L^^) 
[5, 8]. Higher orders are caused by corrections to the leading scaling behavior (for details see [30]). Here we restrict 
ourselves to orders up to 0{L~^). The constant term ao is caused by surface contributions introduced by the considered 
boundary conditions and is absent in the case of periodic and antipcriodic boundary conditions. The numerical data 
and the associated fits are shown in Figure 4 and the results for the fitting parameters are assembled in Table I. 
Interpreting the fit parameter 02 as the Casimir amplitude we estimate the following values: 



Ac.oo = -0.012(1) 
Ac,DO = -0.012(3) 
Ac.oo = -0.012(3) 



(4.7) 



where the numbers in brac;kcts give the error bars. 

As anticipated the numerical values of the Casimir amplitudes lie all together in a small range. They are equal up 
to the third decimal place when we take the adjusted error bar into account. We see that due to the fact that the 
boundary conditions are all in the same surface universality class the evolutions of the excess free energies of the three 
systems are asymptotically equivalent for large L. Looking at the non-universal corrections to the leading scaling 
behavior namely the terms of order 0(1/ L^) and 0(1/ L^) we see that their weights are much higher in the DD and 
DO case compared to the 00 case. 

We checked the accuracy of our method by computing results for the periodic and the antiperiodic case and 
comparing them to the analytical results. Our fit parameters are also shown in Fig. I and the numerical result are 
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compared to the analytic results in Fig. 5. The analytical results for the Casiniir amplitudes [12, 17] are given by 

AcPER = = -0.153... (4.8) 

Cl2(7r/3) C(3) 

Ac,APER = z ? — ^ U.274... . 

3 on 

where CI2 denotes the Clausen function. For the case of periodic boundary conditions we find a perfect agreement 
with the analytic results. For the case of antiperiodic boundary conditions we find an asymptotic agreement for large 
L because there are more corrections due to the inhomogeneity introduced by this boundary condition. 



V. SUMMARY AND CONCLUDING REMARKS 



We showed a new method for the computation of Casimir amplitudes in the large-n limit. This method can be 
applied to the case of slabs with non translationally invariant boundary conditions as well as for those with periodic 
and antiperiodic boundary conditions. We evaluated the Casimir amplitudes for the cases of Dirichlet-Dirichlet, 
open-open and Dirichlet-open boundary conditions and reproduced their asymptotical equivalence in the limit of 
large systems. We find that the non-universal corrections to the leading scaling behavior in the open-open case are 
much smaller than in the other cases. From this we conclude that evaluations of Casimir amplitudes following our 
method in this universality class are simpler in this open-open case. We also confirmed the known analytic result for 
the Casimir amplitudes in the case of periodic and antiperiodic boundary conditions. 
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